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Abstract 

We study the deterministic effects of Raman-induced crosstalk in amplified wavelength division 
multiplexing (WDM) optical fiber transmission lines. We show that the dynamics of pulse ampli- 
tudes in an A^-channel transmission system is described by an A^-dimensional predator-prey model. 
We find the equilibrium states with non-zero amplitudes and prove their stability by obtaining the 
Lyapunov function. The stability is independent of the exact details of the approximation for the 
Raman gain curve. Furthermore, we investigate the impact of cross phase modulation and Ra- 
man self and cross frequency shifts on the dynamics and establish the stability of the equilibrium 
state with respect to these perturbations. Our results provide a quantitative explanation for the 
robustness of differential-phase-shift-keyed WDM transmission against Raman crosstalk effects. 

PACS numbers: 42.65.Dr, 42.81.Dp, 42.81.-1 
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I. INTRODUCTION 



One of the important nonlinear processes affecting pulse propagation in massive wave- 
length division multiplexing (WDAd) optical fiber communication systems is due to inter- 
pulse Raman- induced crosstalk l[ Q]. In this process, which takes place during collisions 
between pulses from different frequency channels (interchannel collisions), energy is trans- 
ferred from high-frequency pulses to low-frequency pulses. It is known that the magnitude of 
the Raman-induced energy exchange in a single interchannel collision is independent of the 
frequency difference between the channels. Consequently, the magnitude of the cumulative 
energy shifts for a given pulse grows with the square of the number of channels, a result that 
is valid for linear transmission js-G], conventional soliton transmission 7 1^, and strongly 
dispersion-managed (DM) soliton transmission ll|]. Therefore, in a 100-channel system, for 
example, the Raman crosstalk effects can be larger by a factor of 2.5 x 10^ compared with 
a two-channel system operating at the same bit rate per channel. 

Early studies of Raman crosstalk in WDM transmission focused on the dependence of the 
energy shifts on the total number of channels jsl , as well as on the impact of energy depletion 
jst and group velocity dispersion I2I, Q on the dynamics. Later on attention turned to 
the combined effects of Raman crosstalk and bit-pattern randomness in the on-off-keying 
(OOK) transmission scheme, and it was found that the probability density function (PDF) 



16| . This finding means that the nth 



of the pulse amplitudes is lognormal 0, 0, 10, \l4 - 
normalized moments of the amplitude grow exponentially with both propagation distance 
and n^. Furthermore, in several studies of conventional soliton transmission it was found 
that the dynamics of the frequency shift is strongly coupled to amplitude djTiamics, and as a 
result, the rath normalized moments of the Raman- induced self and cross frequency shifts also 
grow exponentially with propagation distance and (see Refs. 17H19[|). This intermittent 
dynamic behavior has important practical consequences, by leading to relatively high bit- 
error-rate (BER) values at intermediate and large propagation distances ITHlSj. Raman 
crosstalk effects were also recently investigated in hut-skipped amplified WDM transmission 



2d . |21| , in cable television overlay passive optical networks 22|, |23| , in optical code-division 



multiple-access transmission systems 24|, and in conjunction with four- wave- mixing 25[ |. 

One of the ways to overcome the detrimental effects of Raman crosstalk on massive WDM 
transmission is by replacing the OOK scheme by alternative encoding schemes, which are 
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expected to be less susceptible to these effects. The differential phase shift keying (DPSK) 
scheme, in which all time slots are occupied and the information is encoded in the phase 
difference between adjacent pulses, is one of the promising encoding methods that have 
attracted much interest in recent years H]- DPSK transmission the Raman- induced 
amplitude dynamics becomes (approximately) deterministic, and important questions arise 
regarding the character of this dynamics. One major question concerns the possibility to 
achieve a stable equilibrium state for the amplitudes in all channels. The study reported in 
Ref. jsl demonstrated that this is not possible in unamplified optical fiber lines. However, 
ater experiments showed that the situation is quite different in amplified WDM transmission 
28l . I29I . More specifically, it was found that the introduction of amplification into the 
system significantly reduces the Raman-induced energy shifts. In the present paper we 
suggest a dynamical explanation for this important experimental observation. Moreover, 
we demonstrate the robustness of DPSK transmission against inter-pulse Raman crosstalk 
effects by showing that equilibrium states with non-zero amplitudes in all channels do exist, 
and by proving the stability of the equilibrium states. 

In the present study we consider optical solitons as an example for the pulses carrying 
the information for the following reasons. First, as mentioned above, the Raman- induced 
energy exchange in pulse collisions is similar in linear transmission, conventional soliton 
transmission, and strongly DM soliton transmission. Second, propagation of conventional 
solitons through an optical fiber is described by the nonlinear Schrodinger (NLS) equation, 
which is an integrable model 30| . Due to the integrability of the model, and to the fact that 
optical solitons are stable stationary solutions of the NLS equation, the derivation of the 
model for the Raman-induced amplitude dynamics can be done in a rigorous manner. Third, 
conventional optical solitons have traditionally been considered as excellent candidates for 
information transmission in high-speed optical fiber lines and in all-optical networks 
Furthermore, state-of-the-art transmission experiments already use all-Raman distributed 
amplification 31n37j, which is the most suitable amplification scheme for conventional soli- 
tons. 

Since we consider transmission systems where the pulses in each frequency channel are 
well-separated, intrachannel four-wave-mixing (FWM) effects are negligible. In addition, we 
assume that loss is compensated by distributed Raman amplification. It is a well-known fact 
in soliton theory that in the absence of loss interchannel FWM products completely vanish 
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after the collisions, see Ref. 



. 



301 for theory and Refs. [7|, |38|] for numerical simulations 



with the NLS equation with and without delayed Raman response. We remark that even 
in systems employing lumped amplifiers (and non-overlapping pulses), interchannel FWM 
can be significantly reduced by dispersion-management. This is true both in the soliton 
and in the linear regime 39|. In contrast, the energy shift in a single collision 



regime 
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between two strongly DM solitons is given by an expression with exactly the same form as 



i.e., strong dispersion- management 



the expression obtained for conventional solitons 
does not reduce Raman crosstalk. 

In deriving the model for Raman-induced amplitude dynamics we fully take into account 
pulse walk-off. In addition, we assume that the pulse sequences in all frequency channels 
are deterministic and that the sequences are either infinitely long or are subject to periodic 
temporal boundary conditions. The first setup approximates long-haul transmission, while 
the second one corresponds to a closed fiber loop experiment. We also assume that the 
constant net gain/loss in each channel is determined by the difference between distributed 
amplifier gain and fiber loss. Notice that in this feature our model is fundamentally different 
from the model derived in Ref 5[. Indeed, since in the latter model all channels experience 
net loss, it does not support an equilibrium state with non-zero amplitudes in all channels. 
In contrast, in our model the net gain/loss of some channels is positive while for other 
channels it is negative, which is the underlying reason for the existence of equilibrium states 
with non-zero amplitudes in all channels. 

Our model for the Raman crosstalk dynamics in an A^-channel system consists of a 
system of coupled nonlinear ordinary differential equations (ODEs) for the amplitudes in 
different channels. The system of coupled ODEs can be described in the jargon of population 



dynamics theory as an A^-dimensional predator-prey model [40|]. After obtaining the model 
we look for equilibrium states with non-zero amplitudes in all channels and establish their 
stability with respect to deviations of the initial amplitudes from the equilibrium values. 
We also investigate the dynamic behavior induced by such deviations. In actual optical 
fiber transmission systems pulse dynamics can be influenced by physical processes other 
than Raman crosstalk. It is therefore important to understand the manner in which these 
additional processes perturb the Raman-induced amplitude dynamics. In the present study 
we investigate the effects of three perturbations due to cross phase modulation (XPM), 
Raman self frequency shift (SFS), and Raman cross frequency shift (XFS). For each of 



4 



these physical processes we construct the perturbed model describing amplitude dynamics 
in an A^-channel system and investigate the existence of equilibrium states with non-zero 
amplitudes in all channels. Furthermore, we study the stability of the equilibrium state and 
the typical dynamic behavior of the amplitudes in two-channel systems in the presence of 
each of the three perturbations. 

The rest of the paper is organized as follows. In Section III Al we derive the unperturbed 
model for Raman-induced amplitude dynamics in WDM systems with 2N + 1 channels, and 
in Section III Bl we find the equilibrium states of the model and investigate their stability. 
In Section III CI we validate the predictions of Section IIIBI by numerical simulations. The 
perturbed models with XPM, Raman SFS, and Raman XFS, are studied in Sections IIII Al 
IIIIB| and IIII C| respectively, for two-channel systems. Section [IV] is reserved for conclusions. 
In Appendix |X] we obtain the perturbed models with XPM, Raman SFS, and Raman XFS 
in WDM transmission with 2N + 1 channels. 



II. DYNAMICS OF DETERMINISTIC RAMAN CROSSTALK - UNPERTURBED 
MODEL 



A. Derivation of the model 



Propagation of short pulses of light through an optical fiber in the presence of delayed 
Raman response is described by the following perturbed NLS equation 



1|: 



id,ij + dlij + 2\ij\^ij = -eRijdt\ 



(1) 



where if) is proportional to the envelope of the electric field, z is propagation distance and 
t is time in the retarded reference frame. The term —e^dtlip^^ represents the first order 
approximation for the fiber's delayed Raman response 4l(| and is the Raman coefficient 
42l |. When e^, = 0, the single-soliton solution of Eq. ([T]) in a frequency channel /3 is given 
by 

exp(2X/3) 



ipp{t,z)=ri(3- 



(2) 



cosh(x^) ' 

where Xfs = rip{t-yii- 2(3 z), xp = ap + P{t - yp) + (r^^ - (3^) z, and 77/3, ap and yp are the 
soliton amplitude, phase and position, respectively. 
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Let us describe the main Raman-induced effects on a single collision between a soliton 
from the jth frequency channel and a soliton from the kth frequency channel. We assume 



that eR ^ ~ f^k\ ^ 1, which is 



systems even for adjacent channels 
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;he typical situation in many WDM transmission 



39| . Under this assumption one can show that the 



most important effect of delayed Raman response on the collision is an 0{eR) change in the 



soliton amplitude 



10|,ll5| 



Ar]j = 2eif/(|j - A;|)sgn(/3fc - l3j)vjVk, (3) 



where the coupling constants f{\j — k\) depend on the specific approximation for the Raman 
gain curve. Notice that due to the inclusion of the — k\) factors, Eq. ([3]) does not 
rely on the triangular approximation for the Raman gain curve. If we adopt the triangular 
approximation, we obtain that f{\j — k\) = 1 for any j and k. The effects of the collision in 
order eR/\f3j — f3k\ will be described in Section UTTl where we obtain the perturbed models 
for amplitude dynamics. Since e/j ^ — < 1, effects of order ej^ and higher will 

be neglected. In addition, third order dispersion and self-steepening are neglected since 
the collision-induced effects of these conservative perturbations on the soliton amplitude 
and frequency are of higher order in both the parameter e characterizing the perturbative 



process and (see, e.g., Refs. 
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3)- 



Consider now WDM transmission systems with 2N + 1 channels and frequency difference 
A/3 between adjacent channels. Our model, which takes into account pulse walk-off, is based 
on the following assumptions. (1) The soliton sequences in all channels are deterministic 
in the sense that all time slots are occupied and each soliton is located at the center of a 
time slot of width T. Furthermore, when considering amplitude dynamics, the amplitudes 
are equal for all pulses from the same frequency channel, but are not necessarily equal for 
pulses from different channels. This setup corresponds, for example, to return-to-zero (RZ) 



transmission with differential- phase-shift-keying [46|. (2) The sequences are either (a) in- 
finitely long, or (b) subject to periodic temporal boundary conditions. Notice that setup (a) 
is an approximation for long-haul transmission systems, while setup (b) is an approximation 
for closed fiber-loop experiments. (3) The gain/loss in each channel is determined by the 
difference between distributed amplifier gain and fiber loss. In particular, for some channels 
this difference can be slightly positive, resulting in small net gain, while for other channels 
this difference can be slightly negative, resulting in small net loss. We emphasize that in 
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this feature our model is essentially different from the models that are usually considered in 
studies of Raman amplification schemes, where it is assumed that all pumps experience net 



loss 



15|. 

To obtain the dynamic equation for the amplitude of the jth-channel solitons we note 
that the distance traveled by these solitons while passing two successive solitons in the j — 1 
or j + 1 channels is Azi^^ = T/(2A/3). We denote by zi the location of the Ith collision of 
a given jth-channel soliton with solitons in the j ' + 1 or j — 1 channel. Using Eq. (|3]) and 
summing over all collisions occurring within the interval {zi-i,Zi], where Zi = Zi^i + Az^\ 
we obtain 

N 

r]j{zi_, + Az«) = r]j{zi_i) + g,7]^{zi_,)Azi'^ + 2eK^ {k - j)f{\j - k\)7]^{zi_i)r]k{zi.i).{4) 

k=-N 

The constant gj on the right hand side of Eq. (j4]) is the net gain/loss coefficient for the 
jth channel, which is assumed to be independent of z. Due to the periodicity of the pulse 
sequences the same equation is satisfied by the amplitudes of all jth-channel solitons. More- 
over, the same equations with different j values, where j = — A^, . . . ,N, describe the dynam- 
ics of the soliton amplitudes in all channels. Going to the continuum limit we obtain 

TV 

9j+CY. ik~j)fi\j-k\)vk , (5) 



dr]j 

d7 = 



k=-N 

where C = Ae^Ap/T and j = — A^, . . . ,N. The system ([5]) gives a complete description of 
the dynamics of the r^^'s. Notice that apart from the important fact that some gain/loss 
coefficients can be positive while others can be negative, the system (|5]) is similar to the 
one obtained in Ref. jsj for the Raman-induced amplitude dynamics of continuous waves 
in unamplified WDM transmission. Thus, our model also describes the Raman-induced 
crosstalk dynamics of continuous waves in amplified WDM transmission systems. 



B. Equilibrium states, stability, and conserved quantities 

In optical fiber communication systems it is usually desired to achieve a steady state in 
which the pulse amplitudes in all channels are equal and constant (independent of 2;) 
We therefore look for an equilibrium state of the system ([5]) in the form 77^^'^'' = 77 > for 
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—N < j < N. Setting the right hand sides of ([5]) equal to zero we arrive at 

N 

9, = -Cv Yl ik-j)f{\j-k\). (6) 

k=-N 

Therefore, the gain required to maintain an equihbrium state with equal amplitudes is not 
fiat with respect to frequency. Instead, the net gain/loss coefficient of high-frequency chan- 
nels should be positive, while that of low-frequency channels should be negative. In other 
words, high-frequency channels should be overamplified, whereas low-frequency channels 
should be underamplified compared with the reference (j = 0) channel. Substituting Eq. 

into Eq. (|5]) we arrive at a slightly simpler form of the model, which is convenient for 
analysis and numerical simulations, 

-j^ = Cv, J2(^-3)fi\j-k\){Vk-v)- (7) 

k=-N 

Equation (I7j) can be described in population dynamics terminology as a predator-prey system 
with 2N +1 species 

The equilibrium states of Eq. ([7]) with non-zero amplitudes are determined by 

N 

E - - - ^) = 0, -N<j< N. (8) 

k=-N 

The trivial solution of Eq. ([H]), that is, the solution with rj^f^''^ = rj > for — < k < N, 
corresponds to the equilibrium state of Eq. ([7]) with equal non-zero amplitudes. However, 
due to the anti-symmetry of {k — j)f{\j — k\) with respect to an interchange of j and k, Eq. 
([8]) has infinitely many non-trivial solutions, and these correspond to equilibrium states of 



Eq. ([7]) with unequal non-zero amplitudes 47|]. In the case where the Raman gain curve is 



described by the triangular approximation it is straightforward to show that the non-trivial 
equilibrium states of Eq. ([7]) are determined by the following two equations: 

N N 

5^ = {2N + Yl ^^^^ = 0- (9) 

k=-N k=-N 

Therefore, in this case the equilibrium states lie on a (2A^ — l)-dimensional plane. Thus, for 
a 3-channel system, for example, the equilibrium states lie on the line segment [r], rj, rf) — 
{b - r7)(l/2, -1, 1/2), where < 6 < Sr^. 
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We now prove stability of all equilibrium states of Eq. ([7]) with non-zero amplitudes, 

= 7]^"'^^ > 0,-N < J < N,hy cons 
this purpose we consider the function 



rij = r]^p^ > 0, — < J < by constructing a Lyapunov function for this equation. For 




(10) 



where rj = {rj^j^, . . . ,rij, . . . ,1]]^). Taking the derivative of VL(r}) along trajectories of the 
system (I7j) we obtain 

„ , N _ (eg) < N N 

j=-N j=-N k=-N 

FromEq. © it follows that ^^_^(k-j)fi\j-k\){r^k-v) = Et-Nik- j)fi\j - k\){r^k- 
rj^^i)y Using this relation together with Eq. ( ITTl) and the anti-symmetry of {k — j)f{\j — k\) 
we arrive at 

dV ^ ^ 

j=-N k=-N 

for rjj > 0, — ^ j < N. It is straightforward to show that each term h{'rij) = rjj — rjj^'^^ + 
\n{rij^'^^ /rjj) on the right hand side of Eq. ffTOj) satisfies h{rij) > for any rjj > 0, where 



h{r]j) = only at rij = f^P^ ■ Therefore, ¥^{7)) > for any vector r) for which T]j > for 
~N < j < N, where equality holds only at equilibrium points. Combining this result with 
the result dVi/dz = along trajectories of the system we conclude that V/, is a Lyapunov 
function of the system ([7]) and therefore the equilibrium states rjj = iq^p^ > 0, —N < j < N 



are stable 



49| . We note that since dVi/dz = rather than dVi/dz < this stability 
means that the values of rij{z) are bounded for any z but do not tend to rjj^'^^ for large 
propagation distance z. Thus, the typical dynamics of the amplitudes rij{z) for initial 
conditions that are off the equilibrium point is oscillatory. Notice that the stability is 
independent of the exact specification of the f{\j — k\) values, and therefore the equilibrium 
states are stable irrespective of the specific details of the approximation for the Raman gain 
curve. Notice also that Eq. (ITO!) actually provides K + 1 independent conserved quantities 
for the model ([7]), where K is the dimension of the solution space of Eq. ([8]). For example, 
if the Raman gain is described by the triangular approximation, K = 2N — 1, and there are 
2N conserved quantities in total. 
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C. Numerical solution of the system ([7]) 



In order to check the predictions of the previous subsection we solve Eq. ([7]) numerically 
by employing a fourth-order Runge-Kutta method. As a concrete example for the physical 
parameter values we consider a 101-channel system operating at 40 Gbits/s per channel 
with frequency spacing Az/ = 100 GHz and dimensionless time slot width T = 5. With this 
choice the pulse width is 5 ps, A/3 = tt, eji = 0.0012, and = 50. These parameter values 
are typical for several state-of-the-art massive multichannel transmission experiments, see 



Ref. 39| and references therein. Taking = — Ips^km"^ and 7 = 4W~^km~^ we obtain 
Pq = 10 mW for the soliton peak power. The dimensionless final propagation distance is 
taken as Zf = 200 corresponding to Xf = 10^ km, but the main features of the dynamics 
can already be observed at considerably shorter distances. We also choose 77 = 1, so that 
the trivial equilibrium state is r^j^'^^ = 1, for all j. 

To illustrate the impact of Raman crosstalk on massive WDM transmission we focus 
attention on amplitude dynamics of solitons from faraway channels. Thus, we start by 
considering a two-channel system consisting of the reference channel {j = 0) and the highest- 
frequency channel j = N = 50. Employing Eq. ([7]) to this two-channel system while 
adopting the triangular approximation for the Raman gain curve we obtain 

dV50 



dz 

dm 

dz 



= 50C?75o (1 - 770) , 
SOCr^o (r^so - 1) • (13) 



The solution of Eq. (fT3!) for ?7o > and r^i > is 

r]5oiz) + r]o{z) - In [r]5o{z)] - In [?7o(^)] = (14) 

where the constant k is determined by the values of ?7o(0) and ?75o(0). The z-dependence of 
the amplitudes tiq and 7750 is shown in Fig. [1] (a) for the initial condition 7750(0) = 1.2 and 
^70(0) = 0.9 . It is clearly seen that the amplitudes oscillate about their equilibrium value 
1] = 1. The dimensionless oscillation period is Zp = 41.8, corresponding to 2090 km [50i] . 
Similar oscillatory behavior is observed for other choices of the initial amplitudes. Figure [1] 
(b) shows the phase portrait for the system ( JT3ll . Since all trajectories are closed curves we 
conclude that the dynamics is indeed periodic and that the equilibrium state is stable. 
Next we consider a three-channel system consisting of the channels j = and j = ±50. 
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Within the Raman triangular approximation the amphtude dynamics is described by 



50C?75o (3 - ?7o - 27]- 



dz 



50 ; 



dv-50 



50Cr7_50 (-3 + 27^50 + Vo) • (15) 



dz 

Figure [2] (a) shows the phase portrait. The sohd hne is the hne segment of stable equilibrium 
points (1, 1, 1) — (6— l)(l/2, —1, 1/2), where < 6 < 3. Since all trajectories are closed orbits 
centered about the equilibrium line we conclude that the equilibrium states are indeed stable 
and that the amplitudes exhibit periodic oscillations. This oscillatory dynamics is illustrated 
in Fig. [2](b) for the initial condition 7750(0) = 0.8, 'r]o{0) = 0.9, r/_5o(0) = 1.1. 

As another example for possible dynamic scenarios exhibited by the unperturbed model 
([7]) we study a four-channel system consisting of the channels j = ±16 and j = ±48. As 
we demonstrate below, a four-channel setup represents the simplest case where deviations 
from the triangular approximation for the Raman gain curve lead to the emergence of new 
dynamical features. We consider the following sets of values for the coupling constants 
/(|j — A;|): /(I) = /(2) = 1 and /(3) = 1 + p, where p = for set (1) and p = 0.1 for set 
(2). Thus, set (1) corresponds to the triangular approximation, whereas set (2) represents a 
10% deviation from the triangular approximation for /(3). For these values of the coupling 
constants the dynamics of the amplitude is described by 

dr]i8 



dz 

driw _ 

dz 

dr]-iQ 



32C?748 [(6 ± 3p) - ?7i6 - 2?7_i6 - 3(1 ± p)r]^4s] 
32C?7i6 (2 ± r/48 - r]_iQ - 2r]_^s) , 

= 32C77_i6 (-2 ± 27^48 ± r]iQ - ?7_48) , 



dz 

32C77_48 [-(6 ± 3p) + 3(1 ± p)r]^s + ^Vw + V-w] ■ (16) 



dr]-48 



dz 

We solve the system (1161) numerically with p = and p = 0.1 and with the initial condition 
?748(0) = 1.2, 7/16(0) = 1.1, 77_i6(0) = 0.95, and 77.48(0) = 0.9. The results of our numerical 
simulations are presented in Figs. |3] and |H As can be seen, in both cases the amplitudes 
exhibit oscillatory dynamics. However, in the triangular approximation case {p = 0) there 
is only a single period for the oscillations, whereas in the non-triangular approximation 
case {p = 0.1) two very different oscillation periods are observed. The first period is not 
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significantly different from the oscillation period for p = 0, while the second period is much 
longer than the first one (see Fig. H]). In fact, in order to verify that the amplitude dynamics 
is indeed periodic for p = 0.1 one has to carry out the simulation up to a final propagation 
distance of Zf = 2000. 

Notice that our choice of 40 Gbits/s per channel systems is mainly for reasons of conve- 
nience, since with this value we can illustrate the oscillatory nature of the Raman-induced 
amplitude dynamics already with two and three channels. Since current soliton-based sys- 
tems work at 10 Gbits/s per channel it is useful to examine one example for the amplitude 
dynamics in such systems. We therefore consider a 361-channel system operating at 10 
Gbits/s per channel with frequency spacing Au = 25 GHz and dimensionless time slot 
width T = 5. For these parameter values the pulse width is 20 ps, A/3 = tt, e/j = 3 x 10~^, 
and = 180. Taking /32 = — 2ps^km~^ and 7 = 2W~^km~^ we obtain Pq = 2.5 mW for 
the soliton peak power. The dimensionless final propagation distance is taken as Zf = 25 
corresponding to Xf = 10^ km, but the main dynamical features can already be observed at 
shorter distances. As before we choose 1] = 1, so that the trivial equilibrium state is rjj^''^ = 1, 
for all j. To illustrate the dynamics we consider a seven-channel system consisting of the 
channels j = 0, j = ±60, j = ±120, and j = ±180 and adopt the triangular approximation 
for the Raman gain curve. The corresponding system of equations for amplitude evolution 
is solved numerically with the initial condition 77i8o(0) = 1-2, ?7i2o(0) = 1.05, ?76o(0) = 1-1 
?7o(0) = 1.15, ?7_6o(0) = 0.98, ?7_i2o(0) = 1.1, and ?7_i8o(0) = 0.95. The results are shown 
in Fig. |5l One can see that pulse amplitudes in all channels exhibit oscillatory behavior 
similar to the one observed in Figures [T] and |2] for the 40 Gbits/s per channel system. Thus, 
equilibrium states with non-zero amplitudes in all channels are stable in both 10 and 40 
Gbits/s per channel systems, in agreement with the predictions in Section [11 B[ 

III. DYNAMICS OF DETERMINISTIC RAMAN CROSSTALK - PERTURBED 
MODELS 

We now turn to study the effects of perturbations on the model described by Eq. ([7]). We 
focus attention on the effects of cross phase modulation (XPM), Raman self frequency shift 
(SFS), and Raman cross frequency shift (XFS). It is well-known that each of these processes 
by itself does not change the soliton amplitude. Thus, in the absence of Raman crosstalk the 
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soliton amplitudes are constant and there is no amplitude dynamics to be considered. As we 
show below XPM, Raman SFS, and Raman XFS affect amplitude dynamics by leading to 
2;-dependence of the collision rate. For these reasons, the three processes can be considered 
as perturbations to the model ([7]). For each of the three processes we develop a perturbed 
model for amplitude dynamics and analyze the stability of the equilibrium state and the 
dynamic behavior in the new model. Here we choose to concentrate on a two-channel 
system consisting of channels j = and j = 1 since such treatment is sufficient in order to 
uncover the main changes in the dynamics compared with the unperturbed model. Without 
loss of generality we assume /(I) = 1 for the coupling constant. The general forms of the 
perturbed models in WDM systems with 2N + 1 channels are obtained in Appendix \M 



A. Effects of cross phase modulation 

The XPM-induced position shift experienced by a soliton in the jth channel as a result 



of a single collision with a fcth-channel soliton is given by (see, e.g., Refs. 15|, |38| ) 



_ 4sgn(/3fc - f3j)r]k . . 

Our goal is to obtain a perturbed model for the dynamics of the amplitudes ?7o and r/i in the 
two-channel system, which takes into account the XPM-induced position shift. It is clear 
that the main effect of the position shift ( ITTll is to lead to a change in the inter-collision 
distance ^z^\ that is, Azi^'* is 2;-dependent in the perturbed model. In order to find this 
2-dependence we first write down the equation for the location of the Zth collision of the 
soliton from the zero time slot in the reference channel. Since this collision is with the soliton 
from the — / time slot in channel 1, the collision distance zi is determined by 

yi,-i{zi) = yofl{zi), (18) 

where the first subscript in yi-i stands for the channel and the second subscript represents 
the time slot. Taking into account the difference in group velocities, 2/3i, and the position 
shifts experienced by the solitons in the collision we find 

+ 2/3iA^«(^;_i) - Ar^,{zi_,)ll3l = yo^-i) + (19) 

Solving Eq. ( !T9|) for A^c^^ we obtain 



/^z^^\zi-i) = [2/o,o(^/-i) - yi,-i{zi-i)] + -| + r]^{zi_^)] . (20) 
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From the definition of zi^i and tlie periodicity of the pulse sequences it follows that 
2/o,o(-2«-i) — yi,-i{zi^i) = T. Therefore, Eq. fl20|) can be rewritten as 

T 



213, 



{l + 5[r/o(^/-i) + r/i(2;i-i)]} 



(21) 



where B = 4:/{T(3f). Notice that Azi^\zi^i) > T/(2/3i) due to the fact that the XPM- 
induced position shifts are positive for the reference channel solitons and negative for the 
solitons in channel 1. The equations for the amplitudes rjo and rji at zi = zi^i+Azc^\zi_i) are 
similar in form to Eq. (jl]). Substituting relation (^T^ into the equations for r]o{zi_i + 
and ?7i(z/_i + Azi^^) and going to the continuum limit we obtain 

Cvo 



dr]i 
dz 
drjo 
dz 



91 



9o + 



B ivo + Vi) 

Cm 



(22) 



l + B{r]o + %)_ 

Equation (|22|) represents the perturbed model that takes into account XPM effects. Notice 
that for the typical values T = 5 and /3i = vr, 5 ~ 0.081, i.e., 5^1. Thus, for such values 
the XPM perturbation can be considered as a weak perturbation. 

We look for equilibrium states of Eq. (!22l) that are of the form tjq^'^^ = t]^,''^ = rj > 0. This 
requirement yields the expressions gi = —Qq = Cri/{1 + 2Br]) for the gain/loss coefficients. 
Thus, an important consequence of the XPM-induced position shifts is a change in the values 
of the gain/loss coefficients that are required for maintaining an equilibrium state with equal 
non-zero amplitudes. Taking into account the modified expressions for qq and gi, we can 
rewrite Eq. (122|) as 

V Vo 



-r- = Cvi 

dz 

dz 



1 + 2Br] 
V 



B {riQ + rii) 



J Jl (23) 

l + 2Br] 1 + B{r]o + r]i)\' ^ ' 

To prove stability of the equilibrium point (?7, 77) we look for a Lyapunov function of the 
system ( l23|) in the form 



Vl(?7o,^?i) = r]\n 



+ r/ln 



1 + 2Br] 
B 



In 



l + B{riQ + rii) 
l + 2Bri 



(24) 



Taking the derivative along trajectories of the system fl23|) we arrive at 



dz 



7] drjQ 7] drji 1 + 2Bri d 

Tjo dz Tji dz 1 + B{riQ + rji) dz 



(25) 
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Replacing drjo/dz and dr]i/dz with the right hand sides of Eq. ( I23!l . we find 



-— = -Cr] 
dz 



V , Vi 



2B7] l + B{r]o + r]i) 



-Cr] 



l + 2Br] l + B{r]o + r]i) 



0. (26) 



l + B{r]o + r]i] 

Hence, Vl is constant along trajectories of Eq. fl23|) . In addition, it is straightforward to 
show that Vl attains its minimum at (77, 77). Since Vl^t], rf) = 0, it follows that Vl^tiq, rji) > 
for any (770,771) such that t/q > and 7/1 > 0. Combining these results we conclude tha t Vj, 
is a Lyapunov function for Eq. ( |23l) and that the equilibrium state (77, r/) is a center 48|, |49|. 
Consequently, (7/, 7/) is a stable equilibrium point and deviations of the initial amplitude 
values from 77 lead to oscillatory dynamics of 770(2;) and rii{z). This means that the XPM 
perturbation does not change the stability properties of the equilibrium state. 

To illustrate these conclusions we numerically solve Eq. fl23l) by use of a fourth-order 
Runge-Kutta scheme. For concreteness we consider a 2-channel transmission system oper- 
ating at 160 Gbits/s per channel with time slot width T = 5. In this system the pulse width 
is 1.25 ps and = 0.0048. Notice that WDM transmission at 160 Gbits/s per channel has 
received much attention in recent years both theoretically 5l|] and experimentally 52h54|. 
Taking (32 = — 0.5ps^km^^ and 7 = 4W~^km~^ we obtain Pq = 80 mW for the soliton peak 
power. The dimensionless final propagation distance is taken as Zf = 1280 corresponding 
to propagation over Xf = 8000 km, but the main dynamical features are observed already 
at much shorter distances. Without loss of generality we choose 77 = 1, so that the equilib- 
rium state is (1, 1). Figure [6] shows the 2;-dependence of rji and t/q with the initial condition 
r/i(0) = 1.1 and r/o(0) = 0.9 for A/3 = 2.0 (a) and A/3 = 10.0 (b). The latter values of 
A/3 correspond to interchannel frequency spacing of Au = 250 GHz and Au = 1250 GHz, 
respectively. In both cases the amplitudes exhibit oscillations about the equilibrium value 
in accordance with the prediction of the analytic calculations. Similar oscillatory dynamics 
is observed for other initial conditions and other values of A/3. Furthermore, as can be 
seen from Fig. [3, all solutions curves in the phase plane are closed. Thus, our numerical 
simulations validate the predictions about stability of the equihbrium state in the presence 
of XPM. 
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B. Effects of Raman self frequency shift 



The Raman-induced SFS experienced by the jth-channel sohtons is given by 55l-l57| 



where we exphcitly take into account the couphng of the frequency shift to the amphtude. 
Consider the effects of the Raman SFS on amphtude dynamics in a two-channel system 
consisting of channels and 1. Since in general tiq^z) ^ f]\{z)., the Raman-induced SFS can 
lead to z-dependence of the frequency difference /3io(-2) = — /3o(^), which in turn would 
lead to z-dependence of the inter-collision distance I\z^}^ = T/(2/3io). Therefore, the impact 
of the Raman SFS on amplitude dynamics in a two-channel system can be taken into account 
by replacing the constant frequency difference A/3 by the z-dependent frequency difference 
Pio{z) in the equations for rjQ and r/i, and by using relation (1271) to obtain the dynamic 
equation for /3io. This calculation yields the following system: 

dm ( c 



H - 4) ■ (28) 



dz 15 

Looking for equilibrium states of the system (128|) in the form (77, 77, A/3) we obtain gi = 
—Qq = Crj for the gain/loss coefficients. Comparing this result with the result obtained in 
section HIl for the unperturbed model we see that the Raman-induced SFS does not change 
the values of the gain/loss coefficients that are required for maintaining an equilibrium state 
with equal amplitudes. Using the values gi = —go = Crj we can rewrite the system (l28l) in 
a simpler form, 



d/Sio 



-T^^nivt-Vt). (29) 



dz 15 

We study the stability of the equilibrium state {r],ri,Al3) of (l29l) by linear stability analysis 
and by numerical simulations. Linear stability analysis predicts a bifurcation at A^uf = 
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(16T?7^/15)^/^. For A/3 > A/3;,i/ the Jacobian matrix of fl29p has two purely imaginary 
eigenvalues and one zero eigenvalue, while for A/3 < APtif all three eigenvalues are real. 
In the latter case one eigenvalue is positive, another is negative, and the third one is zero. 
For A/3 = APhif all three eigenvalues are equal to zero. Based on this analysis one might 
suspect that the equilibrium state {ri,ri,Af3) becomes unstable for A/3 < Af3hif. However, 
since (y, n, AB) remains a non-hyperbolic equilibrium point, linear stability calculations 



might fail 48|, |49| and one has to resort to numerical simulations to study stability. 



We perform numerical simulations with Eq. ( l29l) for the same two-channel system that 
was considered in subsection IIIIAI Choosing r/ = 1, the equilibrium state is (1,1, A/3) 
and the bifurcation value predicted by linear stability analysis is Afiuf — 2.31. Figure [8] 
shows the z-dependence of the amplitudes and frequency difference for A/3 = 2.0 (a) and 
A/3 = 5.0 (b). The initial conditions are r/i(0) = 1.001, 170(0) = 0.999, and /3io(0) = 1.999 
in (a), and 771(0) = 1.1, r7o(0) = 0.9, and /3io(0) = 5.1 in (b). One observes that in (a) 
rji{z), rio{z), and (3io{z) tend away from their equilibrium values (1, 1,2), while in (b) ?7i(z), 
770(2;), and (3io{z) oscillate about their equilibrium values (1, 1, 5). Similar dynamic behavior 
is obtained for other initial conditions in the neighborhood of (1,1,2) or (1,1,5). Based 
on these observations we conclude that the equilibrium state (1,1,2) is unstable, whereas 
(1,1,5) is stable, in agreement with linear stability analysis. This conclusion is further 
supported by the corresponding phase portraits that are shown in Fig. |9l For the system 
with A/3 = 2.0 the trajectories starting in the vicinity of (1, 1,2) tend away from (1, 1,2). 
In contrast, for the system with A/3 = 5.0 trajectories in the vicinity of (1, 1,5) are closed 
curves centered about (1, 1,5). 

To further investigate the stability properties of a generic equilibrium point (7/, 7/, A/3) 
we perform detailed analysis of numerical simulations data for t] and A/3 values close to the 
bifurcation line predicted by linear stability calculations: A^bif ~ 2.311]^'^. Figure ITO] is the 
bifurcation diagram obtained by this analysis. It is seen that the bifurcation line obtained 
by numerical solution of Eq. ( 129|) closely agrees with the line predicted by linear stability 
computations despite of the fact that (77, rj, A/3) is a non-hyperbolic equilibrium state. Notice 
that the points in the rj — A/3 plane that are above the bifurcation line are stable, while 
those that are below the bifurcation line are unstable. Thus, transmission in the two-channel 
system described in subsection IIIIAI becomes unstable for frequency spacing values smaller 
than 2.317/1-^ 
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C. Effects of Raman cross frequency shift 



The second effect of delayed Raman response on a single collision between a fth-channel 



soliton and a fcth-channel soliton is an OjeulMij — Pik\) frequency shift |7l. lOHll 
Raman- induced XFS is given by 1^, 17, 18| 
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This 



We now consider the impact of this frequency shift on amplitude djTiamics in a two-channel 
system. Since in general the amplitudes t^q and r^i vary with ^, the frequency shifts might 
lead to z-dependence of the frequency difference /3io. In order to construct the perturbed 
model describing the effects of the Raman XFS we first obtain an equation for the dynamics 
of /3io. Using Eq. fl5Ul) we find dPo/dz = — 16e/j?7gr7i/(3T) and df3i/dz = — 16e/jr7o^i/(3T). 
Combining these relations with the definition of /3io we arrive at 



d(3io 
dz 



-IQenVoViiVi-Vo) /i^T). (31) 



In addition, we replace the constant frequency difference A/3 by the z-dependent frequency 
difference /3io{z) in the equations for rjo and rji to obtain 

'^^^ Vo (go + -^PioVi ) • (32) 



dz V 

Equations (13T|) and ( 132|) represent the perturbed model for amplitude dynamics in the pres- 
ence of the Raman XFS. 

We look for equilibrium states of the model in the form [i], t], A/3) and find gi = —go = Crj 
for the gain/loss coefficients. Therefore, the Raman XFS does not modify the values of the 
gain/loss coefficients required for maintaining an equilibrium state with equal amplitudes. 
Substituting gi = —go = Crj into Eq. ( l32l) we obtain the following simpler form of the 
model: 

^ = Cr/i {rj - /3ior/o/A/3) , 

^ = Cr,o {-V + PioVi/m , 
dz 

= -IQtRijorii ivi - Vo) I (3r). (33) 
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It is possible to show that Vjs:(?7i, r^o, /Sio) = Sr^o^^i/S + /3fo is a conserved quantity for the 
system f p3|) . However, this does not guarantee the stabihty of the equihbrium point. Linear 
stabihty analysis predicts bifurcation at Af3bif = (S/SY^'^r]. For A/3 > A/J^jj the Jacobian 
matrix of fl55]) has two purely imaginary eigenvalues and one zero eigenvalue, while for 
A/3 < A/3f,jj all three eigenvalues are real. In the latter case one eigenvalue is positive, 
another is negative, and the third one is zero. For A/3 = APbif all three eigenvalues are 
equal to zero. Based on this analysis one suspects that the equilibrium state [t], t], A/3) 
becomes unstable for A/3 < Af3bif. 

Since (77, 77, A/3) is a non-hyperbolic equilibrium point linear stability calculations might 



lead to erroneous conclusions 
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49|. We therefore investigate the stability of {rj, r], A/3) by 



numerical simulations with Eq. ([33]) • For concreteness we consider the two-channel system 
described in subsection IIII A[ Choosing r] = 1, the equilibrium state is (1,1, A/3) and the 
bifurcation value predicted by linear stability analysis is Af3bif — 1.63. Figure [TT] shows 
the dynamics of the amplitudes and frequency difference for A/3 = 1.5 (a) and A/3 = 5.0 
(b). The initial conditions are r]i{0) = 1.001, r]o{0) = 0.999, and /3io(0) = 1.499 in (a), and 
?7i(0) = 1.2, 770(0) = 0.9, and /3io(0) = 5.05 in (b). It is observed that in (a) rii{z), riQ{z), and 
/3iq{z) tend away from their equilibrium values (1, 1, 1.5), while in (b) rii{z), tiq^z), and f3io{z) 
oscillate about the equilibrium values (1, 1,5). Additional numerical simulations with other 
initial conditions in the vicinity of (1, 1, 1.5) or (1, 1, 5) show the same dynamical behavior. 
We therefore conclude that the equilibrium state (1,1,1.5) is unstable, while (1,1,5) is 
stable, in agreement with linear stability analysis. This conclusion is also supported by the 
corresponding phase portraits. 

To check the stability of a generic equilibrium state (77, 77, A/3) we carefully analyze re- 
sults of numerical simulations with Eq. fl33|) for different rj and A/3 values. We pay special 
attention to the region in the 7/- A/3 plane that is in the close neighborhood of the bifurcation 
line A/3f,j/ ~ 1.631], predicted by linear stability computations. Figure [12] shows the bifurca- 
tion diagram that is obtained by our analysis. It is seen that the bifurcation line obtained 
by numerical simulations is in good agreement with the line predicted by linear stability 
calculations, despite of the non-hyperbolic character of (77,77, A/3). Notice that equilibrium 
points lying below the line AP = 1.637] are unstable, while those lying above it are stable. 
As a practical consequence we note that transmission in the two-channel system described 
in subsection IIII Al is unstable for frequency spacing values smaller than 1.637]. 
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IV. CONCLUSIONS 



We studied the deterministic effects of inter-pulse Raman crosstalk in amplified WDM 
optical fiber transmission systems. We considered conventional optical solitons as an ex- 
ample for the pulses carrying the information and assumed that the pulse sequences in 
all frequency channels are deterministic and that the sequences are either infinitely long 
or arc subject to periodic temporal boundary conditions. The first setup approximates 
return-to-zero (RZ) differential-phase-shift-keyed (DPSK) long-haul transmission, while the 
second one corresponds to RZ DPSK closed fiber loop experiments. We assumed in addi- 
tion that the constant gain/loss in each frequency channel is determined by the difference 
between distributed amplifier gain and fiber loss. Under these assumptions we showed that 
the dynamics of pulse amplitudes in an A'^-channel transmission line is described by a sys- 
tem of coupled nonlinear ordinary differential equations (ODEs), having the form of an 
A'"-dimensional predator-prey model. We calculated the gain/loss coefficients required for 
maintaining an equilibrium state with equal non-zero amplitudes in all channels, and showed 
that high-frequency channels should be over amplified, while low-frequency channels should 
be underamplified compared with the middle (reference) channel. This means that the net 
gain/loss profile should not be fiat with respect to the frequency. With these values of the 
gain/loss coefficients we proved stability of equilibrium states with non-zero amplitudes in 
all channels by constructing Lyapunov functions for the system of ODEs. The stability 
was found to be independent of the exact details of the approximation for the Raman gain 
curve. Furthermore, since the Lyapunov functions are conserved quantities for the system, 
typical dynamics of the amplitudes for initial conditions that are off the equilibrium points 
is oscillatory. 

In an actual optical fiber line Raman crosstalk is not the only process impacting pulse 
dynamics. It is therefore important to understand the manner in which other physical pro- 
cesses perturb the Raman-induced amplitude dynamics described above. In this study we 
concentrated on the effects of three perturbations due to cross phase modulation (XPM), 
Raman self frequency shift (SFS) and Raman cross frequency shift (XFS). For each of these 
physical processes we constructed the corresponding perturbed model for an A^-channel sys- 
tem and studied the dynamics in a two-channel system. For XPM-pcrturbcd two-channel 
transmission we found that the gain/loss coefficients required for maintaining an equilib- 
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rium state with equal non-zero amplitudes are smaller compared with the values in the 
unperturbed case. This is explained by noting that the XPM-induced position shifts tend 
to increase the inter-collision distance and thus to decrease the rate of collisions. In con- 
trast, the stability of the equilibrium state with equal non-zero amplitudes is not changed 
by XPM, i.e., the equilibrium state remains a center. For two-channel systems perturbed by 
Raman SFS or Raman XFS we found that the values of the gain/loss coefficients required 
for maintaining the equilibrium state are the same as in the unperturbed case. However, 
the stability properties of the equilibrium state change as the frequency difference between 
the channels is decreased or increased, i.e., the system undergoes a bifurcation. The bi- 
furcation curves are given by Af3bif = (16T?7^/15)^/^ for the perturbed model with Raman 
SFS and A(3iyif = (8/3)^/^?7 for the perturbed model with Raman XFS. In both models, for 
a fixed value of 77, two-channel transmission with A/3 > ^Puf is stable, while two-channel 
transmission with A/3 < Al3f,if is unstable. We therefore conclude that the Raman-induced 
interplay between amplitude dynamics and frequency dynamics sets a bound on the smallest 
frequency spacing for stable transmission. 

In summary, our study provides a quantitative explanation to the stability of WDM 
DPSK transmission against Raman crosstalk effects. This stability was demonstrated in 
experiments in a closed fiber loop [23!]. The stable behavior of Raman-induced amplitude 
dynamics in DPSK transmission is very different from the intermittent dynamic behavior 
exhibited by pulse pa.a.ete. . o,.off-keyed ,OOKLt_.on due to the iutetplay 



between Raman crosstalk and bit-pattern randomness [4. l6l. I17H19[|. This different dynamic 



behavior is an important advantage of DPSK transmission over OOK transmission. 



Appendix A: Perturbed models in WDM transmission with 2N + 1 channels 

In this appendix we derive the perturbed models for Raman-induced amplitude dynamics 
in the presence of XPM and Raman SFS and XFS for WDM transmission lines with 2N + 1 
channels. 
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1. Cross phase modulation 



The rate of collisions of a soliton from the jth channel with solitons from the the fcth 
channel in the unperturbed model is 

R], = \k-j\/Azi'\ (Al) 

where Azi^'^ = T/(2A/3) is constant. Thus, Eq. (j3]) for the change of the soliton amplitude 
in the interval zi^i + A2;c"'^^] in the unperturbed model can be written as 

7]^{zi_i + A^i)) = 7],{zi.^) + gj7]j{zi_^)Az^^^ 

N 

+2enAzi''>J2 R]k^gn{k - j) f {\j - k\)vAzi^r)vk{zi-i). (A2) 

k=-N 

In the perturbed transmission system, Eq. flAip is replaced by 

RU^) = \k-j\/^zl^l{z), (A3) 

where the 2r-dependent inter-collision distance Az^j^ is affected by the XPM-induced position 
shifts. Replacing K^f^ by in Eq. (]A2p while employing relation (1A3P we arrive at 

r],{zi^i + A^i)) = 7]j{zi^,) + g,r],{zi^{)Az^p 



+2enAz^}^Y^ f{\j - k\)v,{zi.,Mzi-i). (A4) 

k=-N ^^cjki^l-l) 

In order to find an expression for Az^j|,(z) we write down an equation for the location of the 
collision of the soliton from the zeroth time slot in the jth channel with the soliton from the 
l{j — k) time slot in the fcth channel: yk,i{j-k){zi) = yj^o{zi). Taking into account the different 
group velocities and summing over all XPM-induced position shifts during the collisions we 
arrive at the following generalization of Eq. (fT9l) : 

4 , \m - k\sgn{f3m - f3k){l - Smk) 



yk,iu-k){zi-i) + 2kAf3Az^^^{zi_i) + -^^^ J2 Vm{zi-i)- [m-kf 



where 6ij is the Kronecker delta function. Solution of Eq. flASp for A2;^j|.(2;;_i) yields 
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^mk 1 ^mj 



m — k m — j 



(A5) 



(A6) 



Taking the continuum limit in Eqs. (lA4p and (1A6P we obtain 



drjj 
dz 



Vj 



N 



k 



(AT) 



Equation (IA7p together with Eq. flA6p (with zi^i replaced by z) represent the perturbed 
model for Raman-induced amplitude dynamics in the presence of XPM in WDM transmis- 
sion systems with 2N + 1 channels. 



2. Raman self and cross frequency shifts 

Consider the perturbed model with the Raman- induced SFS. The change in the amplitude 
of a jth-channel soliton within the interval {zi^i,zi^i + Azc^^] is still given by Eq. 
However, now the z-dependent inter-collision distance Az^^jl is given by 



= - j)T/mjM), (A8) 
where Pkj{z) = (3k{z) — /3j{z). The dynamics of (3kj is governed by 

Substituting relation (lASP into Eq. flA4p and going to the continuum limit we obtain 

TV 

+ fi\3-k\)Pk,Vk . (AlO) 

k=-N 

Equations (1A9P and (lAlOp with — ^ j,k < N describe the Raman-induced amplitude 
dynamics in the presence of Raman SFS in transmission systems with 2N + 1 channels. It is 
straightforward to show that the gain/loss coefficients gj that are required for maintaining 
an equilibrium state with equal non-zero amplitudes in all channels are given by Eq. (|6]), 
that is, the values of these coefficients are not modified by the Raman SFS. 

Turning to the perturbed model with Raman XFS we observe that amplitude dynamics 
is described by Eq. flAlOp . To obtain the dynamic equation for the frequency difference Pkj 
we first compute the change in the frequency experienced by a soliton in the jth channel 
within the interval {zi-i,zi-i + Azc^^]. Employing Eq. ( 130|) and summing over all collisions 
we arrive at 



dr]j 

d7 = ^^- 



+ Az«) = - ^V]izi-i) E Vkizl.^)il - 6,,). (All) 

^ k=-N 
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The continuum limit of Eq. (lAlip is 



dl3j 



3T 



N 



dz 



(A12) 



k=-N 



and therefore the dynamics of Pkj{z) is governed by 



dl3kj 
dz 



16eR 

3r 



N 



VkVj{Vk-Vj) + Yl ^"^^^k 



v!)i^-5mk)il-5mj) . (A13) 



m=-N 



Thus, the perturbed model for Raman-induced amplitude dynamics in the presence of Ra- 
man XFS is given by Eqs. ( lAlOl) and ( 1A13|) . where — < j,k < N. A simple calculation 
shows that the values of the gain/loss coefficients required to maintain an equilibrium state 
with equal non-zero amplitudes in all frequency channels are the same as in the unperturbed 
model. 
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List of Figure Captions 

Fig. 1. Amplitude dynamics for the two-channel system described by Eq. (fT3|) . (a) The z- 
dependence of 1750 (solid curve) and r^o (dotted curve) with the initial condition ^750(0) = 1.2 
and ?7o(0) = 0.9. (b) The corresponding phase portrait. 

Fig. 2. Amplitude dynamics for the three-channel system described by Eq. (fT5|) . (a) The 
phase portrait. (1,1)1) is the second equilibrium point from the right, (b) The 2;-dependence 
of r^so (solid curve), r^o (dotted curve), and r^.so (dashed curve), with the initial condition 
7750(0) = 0.8, r7o(0) = 0.9, r7_5o(0) = 1.1. 

Fig. 3. The z-dependence of pulse amplitudes for the four-channel system described by Eq. 
ffTB]) for p = (a), and p = 0.1 (b). The initial condition is 7743(0) = 1.2, ?7i6(0) = 1.1, 
?7_i6(0) = 0.95, and 77_48(0) = 0.9. The solid, dashed, dashed-dotted, and dotted curves 
represent 7743(2;), r]x%{z)-, r]_x%{z\ and r]_^^{z)^ respectively. 

Fig. 4. The 2;-dependence of the amplitudes of pulses in the j = 48 channel for p = (a), 
and p = 0.1 (b). The initial condition is the same as in Fig [3] and the final propagation 
distance is = 2000. 

Fig. 5. The z-dependence of pulse amplitudes for a seven-channel system operating at 10 
Gbits/s per channel. The initial amplitudes are 77130(0) = 1.2, 7/120(0) = 1.05, 7/6o(0) = 1.1 
r/o(0) = 1.15, r/_6o(0) = 0.98, r/_i2o(0) = 1.1, and r/_i3o(0) = 0.95. The solid, dashed, dotted, 
and dashed-dotted curves in (a) represent 77130(2;), 7759 (z), 77_eo(-2), and 77_i8o(2;), respectively. 
The solid, dashed, and dotted curves in (b) correspond to 77120(2;), 770(2;), and 77. 120 (-2)- 
Fig. 6. The 2;-dependence of pulse amplitudes for the XPM-perturbed two-channel system 
described by Eq. ^ with A/3 = 2.0 (a) and A/3 = 10.0 (b). The initial condition 
is r7i(0) = 1.1 and 770(0) = 0.9. The solid and dotted lines represent ?7i(2;), and ?7o(^), 
respectively. 

Fig. 7. The phase portrait for the XPM-perturbed two-channel system described by Eq. 
(1231) with A/3 = 2.0 (a) and A/3 = 10.0 (b). 

Fig. 8. The 2;-dependence of soliton amplitudes and frequency difference for a two-channel 
system perturbed by the Raman SFS [Eq. (12^ ]. (a) The dynamics with A/3 = 2.0 and 
initial condition 771(0) = 1.001, 770(0) = 0.999, and /3io(0) = 1.999. (b) The dynamics with 
A/3 = 5.0 and initial condition ?7i(0) = 1.1, 770(0) = 0.9, /3io(0) = 5.1. The sohd, dashed and 
dotted lines represent ?7i(z), 77o(-2), and /3io(2;), respectively. 
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Fig. 9. The phase portraits for the two-channel system perturbed by the Raman SFS with 
A/3 = 2.0 (a) and A/3 = 5.0 (b). 

Fig. 10. The bifurcation diagram for the two-channel system perturbed by the Raman SFS. 
The squares correspond to the bifurcation values obtained by numerical solution of Eq. fl2^ . 
while the solid line is a fit of the form A/3fejj = 2.30?7^'^^ for the numerical data. The dotted 
line represents the prediction of linear stability analysis: Afiuf — 2.31?7^'^. 
Fig. 11. The 2;-dependence of soliton amplitudes and frequency difference for a two-channel 
system perturbed by the Raman XFS [Eq. (!33|) ]. (a) The dynamics with A/3 = 1.5 and 
initial condition r]i{Q) = 1.001, r]o{0) = 0.999, and /3io(0) = 1.499. (b) The dynamics with 
A/3 = 5.0 and initial condition r7i(0) = 1.2, 770(0) = 0.9, /3io(0) = 5.05. The solid, dashed 
and dotted lines represent ?7i(z), rjo^z), and f3io{z), respectively. 

Fig. 12. The bifurcation diagram for the two-channel system perturbed by the Raman XFS. 
The squares correspond to the bifurcation values obtained by numerical solution of Eq. (!33|l . 
while the solid line is a fit of the form AfSuf = 0.04 + l.GAi] for the numerical data. The 
dotted line is the prediction of linear stability analysis: A(3bif — l.GSr/. 
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(a) (b) 




FIG. 1: Amplitude dynamics for the two-channel system described by Eq. (jl3p . (a) The z- 
dependence of 7/50 (solid curve) and r/o (dotted curve) with the initial condition ??5o(0) = 1.2 and 
r/o(0) = 0.9. (b) The corresponding phase portrait. 

(a) (b) 
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FIG. 2: Amplitude dynamics for the three-channel system described by Eq. (jlSh . (a) The phase 
portrait. (1,1,1) is the second equilibrium point from the right, (b) The z-dependence of r/50 
(solid curve), ??o (dotted curve), and 77-50 (dashed curve), with the initial condition r/5o(0) = 0.8, 
r?o(0) = 0.9, r?_5o(0) = 1.1. 




FIG. 3: The z-dependence of pulse amplitudes for the four-channel system described by Eq. (fT6ll 
for p = (a), and p = 0.1 (b). The initial condition is r?48(0) = 1.2, r?i6(0) = 1.1, r/-i6(0) = 0.95, 
and r/-48(0) = 0.9. The solid, dashed, dashed-dotted, and dotted curves represent ??48(-z), ??i6(-2), 
r\-\%{z)^ and ry_48(2;), respectively. 
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FIG. 4: The z-dependence of the ampHtudes of pulses in the j = 48 channel for p = (a), and 
p = 0.1 (b). The initial condition is the same as in Fig [3] and the final propagation distance is 
Zf = 2000. 
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FIG. 5: The z-dependence of pulse amplitudes for a seven-channel system operating at 10 Gbits/s 
per channel. The initial amplitudes are ??i8o(0) = 1-2, r/i2o(0) = 1.05, r/6o(0) = 1.1 %(0) = 1.15, 
^-60 (0) = 0.98, ??-i2o(0) = 1.1, and r7_i8o(0) = 0.95. The solid, dashed, dotted, and dashed-dotted 
curves in (a) represent ?/i8o(-2), Veoiz), r]-Qo{z), and r]^isQ{z), respectively. The solid, dashed, and 
dotted curves in (b) correspond to r]i2o{z), r]Q{z), and r/_i2o(-z)- 
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FIG. 6: The z-dependence of pulse amplitudes for the XPM-perturbed two-channel system de- 
scribed by Eq. ([23]) with A/3 = 2.0 (a) and A/3 = 10.0 (b). The initial condition is ryi(O) = 1.1 and 
r/o(0) = 0.9. The solid and dotted lines represent r]i{z), and r]o{z), respectively. 
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FIG. 7: The phase portrait for the XPM-perturbed two-channel system described by Eq. (j23p with 
A/3 = 2.0 (a) and A/3 = 10.0 (b). 




FIG. 8: The z-dependence of sohton amphtudes and frequency difference for a two-channel system 
perturbed by the Raman SFS [Eq. (j29p ]. (a) The dynamics with A/3 = 2.0 and initial condition 
r/i(0) = 1.001, T/o(0) = 0.999, and /3io(0) = 1.999. (b) The dynamics with A/3 = 5.0 and initial 
condition 771(0) = 1.1, ?/o(0) = 0.9, /3io(0) = 5.1. The solid, dashed and dotted lines represent 
i]o{^)j ^-iid f3io{z), respectively. 
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FIG. 10: The bifurcation diagram for the two-channel system perturbed by the Raman SFS. The 
squares correspond to the bifurcation values obtained by numerical solution of Eq. (I29p . while the 
solid line is a fit of the form A/J^jj = 2.30r/^'^^ for the numerical data. The dotted line represents 
the prediction of linear stability analysis: A/35jj ~ 2.31r/^'^. 




FIG. 11: The ^-dependence of soliton amplitudes and frequency difference for a two-channel system 
perturbed by the Raman XFS [Eq. (j33p ]. (a) The dynamics with A/3 = 1.5 and initial condition 
r/i(0) = 1.001, ??o(0) = 0.999, and /3io(0) = 1.499. (b) The dynamics with A^ = 5.0 and initial 
condition 771(0) = 1.2, r/o(0) = 0.9, /3io(0) = 5.05. The solid, dashed and dotted lines represent 
?7i(z), r]o{z), and /3io(z), respectively. 
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FIG. 12: The bifurcation diagram for the two-channel system perturbed by the Raman XFS. The 
squares correspond to the bifurcation values obtained by numerical solution of Eq. (I33p . while the 
solid line is a fit of the form A/J^jj = 0.04 + 1.64// for the numerical data. The dotted line is the 
prediction of linear stability analysis: A/3fcjj ~ 1.6877. 
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